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CHAPTER  I 

IHTRCDUCTION 

Purpose  Of  The  Study: 

Rheology  is  the  science  of  the  deformation  and  flow  of  matter. 
Rheology  today  is  at  a  turning  point.  Many  of  the  problems  of  the  past 
have  received  satisfactory' answers .  The  general  phenomenological  theory 
for  time- independent  flow  of  all  types  in  capillary  and  rotating-cylinder 
viscometers  has  been  developed. 

In  spite  of  the  fact  that  many  advancements  have  been  made  in  the 
field  of  rheological  measurements;  there  are  very  few  viscometers 
commercially  available  which  can  determine  the  viscosity  of  fluids  ranging 
up  to  10   poises  which  will  produce  satisfactory  results. 

In  view  of  this,  it  was  proposed  to  design  a  viscometer  for  determining 

10 
viscosity  up  to  10   poises.  On  the  basis  of  theoretical  analysis  of  Goudy 

14 

and  Kirmser,   a  viscometer  was  designed  and  fabricated  at  Kansas  State 

University.  The  author  used  this  viscometer  for  purpose  of  determining 
viscosity  of  various  fluids.  The  purpose  of  this  report  is  to  present  the 
details  of  the  experiment  performed  and  data  obtained  thereby.  The  data 
was  analysed  statistically  using  computer  program  to  obtain  the  results. 
In  order  to  compare  these  results  the  author  used  Brookfield  Viscometer 
also.   The  details  of  results  of  both  has  also  been  discussed. 


Review  Of  Literature: 

What  is  Viscosity: 

Viscosity  is  the  property  of  a  fluid  (  liquid  or  gas  )  that  mainly 
characterizes  its  flow  behaviour  and  by  virtue  of  which  it  offers  resistance 
to  shear  stress.   Viscosity  defines  the  internal  friction  betv?een  the 
molecules  of  a  fluid.  Whenever  any  layer  of  a  fluid  is  caused  to  move  over 
another,  a  frictional  resistance  is  offered  by  the  accompanying  layers  and 
tend  to  be  carried  along  too.  Due  to  this  resistance  are  developed 
velocity  differences  within  the  fluid  and  forms  the  basis  of  the 
quantitative  measure  of  viscosity. 

The  measurement  of  viscosity  is  of  considerable  importance  in  both 
industrial  production  and  fundamental  science.  Viscosity  helps  in 
determining  the  forces  to  be  overcome  when  fluids  are  used  In  pipelines, 
bearings  etc.   It  also  helps  in  controlling  the  flow  of  liquid  in  processes 
like  spraying,  injection  moulding,  extrusion  and  surface  coating.   It  has 
an  important  bearing  with  mixing  and  heat  transfer  characteristics  of 
fluids.   Viscosity  influences  the  size,  shape  and  arrangement  of  the 
molecules  in  fluids  and  as  such  it  is  of  great  importance  to  a  chemist. 

Newton's  Law  of  Viscosity: 

It  states  that  for  a  given  rate  of  angular  deformation  of  fluid  the 
shear  stress  is  directly  proportional  to  the  viscosity. 
Mathematically 

shearing  stress 
Viscosity  (U)   « __~„ 


rate  of  shearing  strain 


Consider  two  parallel  plates  of  surface  area  A  filled  with  a  fluid  as 

Shown  in  Fig.  1.   A  tangenlial  force  F  is  applied  to  the  top  plate  to  nove 

it  at  a  constant  velocity  U  parallel  to  the  stationary  botto;i  plate. 

V 
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Stationary  Plate 


Fig.  1.  Flow  between  parallel  plates 
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Neglecting  gravitational  forces 


A     b 


or 

I  -Ml  -  4 
A     b 

Where  the  constant  of  proportionality  is  defined  as  the  coefficient  of 

viscosity  (/£). 


If  the  velocity  distribution  is  as  shown  in  Fig.  2  then  UA 


Therefore 


du 

3y 


•/  ~M 


du 

dy 


Fig.  2.  Velocity  distribution 


Fluids  which  obey  Newton's  Law  of  viscosity  are  classified  as  Newtonian, 

In  Newtonian  fluid  there  is  a  linear  relation  between  the  magnitude  of 
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applied  shear  stress  and  the  resulting  rate  of  de Cor nation  as  shown  in  the 

i 
Fig.  3, 
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Fig.  3.  Rheological  diagram 


while  the  non-IIewtonian  fluid  follows  a  non- linear  path, 
discussion  is  United  to  Newtonian  fluid  only. 


Primarily  this 


Measurement  of  Viscosity: 

Experimently  the  viscosity  of  a  fluid  can  be  :. -  ith  an 

Instrument  called  viscometer.      There  are  a  variety  i  _ters   available 

and  can  be  classified   into  three  main  types: 

1.  Capillary  Viscometer 

2.  Rotational  Viscometer  * 

3.  Falling  Ball  Viscometer 

Discussion  of  each  type  in  detail,  follows. 

Capillary  Viscometers: 
Introduction: 

Probably  the  first  scientific  experiment  in  which  -    miliary  or  tube 
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was  used  to  measure  flow  was  made  by  Hagen  in  1839.     i  ,^v.   shortly  then 

after  by  Poiseuille's  work.  Poiseuille  studied  cap11'-'  few  problems  in 

order  to  understand  better  the  circulation  of  blood  :  < i  capillary 

vessels  in  the  human  body.  He  discovered  the  relate*  ;  liown  as  the  Hagen- 

Poiseulle  Law  )  between  flow  rate  and  pressure  drop  ton  I  J. llary  flow 

after  experimenting  fortunately,  with  water  instead  r<£  "  ■■!.      Probably  he 

would  not  have  arrived  at  this  result  had  he  used  b'-ro  .  non-Newtonian 

fluid!   This  discovery  was  the  foundation  of  capilL-r_  ..c  ometery. 
Following  Poiseuille,  Wiederman  and  later  Hagenbach  <ie '   -1  a  theoretical 

formula  for  Poiseuille's  discovery  based  on  Newton'c  uion  of  viscosity. 

Until  1890,  when  Couette  devised  a  new  method  based  on  "tern  of  two 

concentric  cylinders,  capillary  flow  was  the  only  vi  ed  technique  in 
viscometry. 

Capillary  viscometers  are  the  most  commonly  use  for  measurements 

on  Newtonian  liquids.  They  are  comparatively  simple  lexpensive,  they 


require  only  a  snail  quantity  of  Lest  liquid,  temperature  control  is  easy 
and  a  full  mathematical  treatment  is  possible.   In  general,  the  liquid  is 
made  to  flow  through  a  capillary  tube  under  a  known  pressure  difference 
and  the  rate  of  flow  is  measured,  usually  by  noting  the  time  taken  for  a 
given  volume  of  the  liquid  to  pass  a  graduation  mark.   In  certain  types  of 
instrument  the  liquid  is  forced  through  the  capillary  at  a  predetermined 
rate  and  the  pressure  drop  thereby  produced  across  the  capillary  is 
measured.  A  complete  capillary  viscometer  consists  of  five  essential  parts: 

1.  Fluid  reservior 

2.  Capillary  of  known  dimensions 

3.  A  unit  for  controlling  and  measuring  the  applied  pressure 
(i.     A  unit  for  determining  flow  rate,  and 

5.  A  unit  for  controlling  temperature 

Commercially  available  capillary  viscometers  may  be  divided  into  three 
main  types. 

a.  The  cylinder-piston  (  or  plunger  )  variety 

b.  Glass  capillary  units,  and 

c.  Orifice  viscometers 

a.  Flow  in  a  Capillary  Tube  (  Poiseuille's  Equation  ) 

Consider  a  cylinderical  capillary  of  radius  'a»  and  length  •  1*  with  a 
pressure  difference  P  between  the  ends.   It  is  assumed:  (i)  that  the  flow 
is  every  where  parallel  to  the  axis  of  the  tube;  (ii)  that  the  flow  is 
steady,  i.e.  there  is  no  acceleration  of  the  liquid  at  any  point;  (iii) 
that  there  is  no  slip  at  the  walls,  i.e.  the  liquid  in  contact  with  the 
walls  of  the  capillary  is  stationary;  (iv)  that  the  liquid  is  Kewtonian, 
i.e.  the  ratio  of  the  shearing  stress  F  to  the  rate  of  shear  D  is  a  constant^ 


the  viscosity  coefficient. 

Lot  tho  velocity  of  the  liquid  at  a  distance  r  from  the  axis  of  the 
capillary  be  v 


Then 


dv 
D  ■  velocity  gradient  =  -  — 

dr 


< 
Fig.  4.  Flow  in  a  Capillary  Tube 


Consider  the  forces  acting  on  a  cylinderical  element  of  length  1  and 
co-a:cial  with  the  capillary.  That  due  to  the  applied  pressure  P  is  #r2P 
acting  in  the  direction  of  notion.  That  due  to  the  viscous  resistance  of 
the  surrounding  liquid  is  F  x  area  -  F  x  2^rl.  Since  there  is  no 
acceleration  of  the  liquid,  these  forces  must  balance,  so  that 

F  x  2  j/rl  =  jfr2? 
or 

F-* 

21 


Newtonian  behaviour  requires  that 

F  »^D,  so  that 

£v  „  .  Pr  ■ 
dr     2^1 

Thus  the  shearing  stress  and  rate  of  shear  are  directly  proportional 
to  r  and  are  zero  on  the  axis  of  the  capillary.  • 

Integrating  the  above  equation,  using  the  condition  that  at  the  wall 
of  the  capillary, 

r  ■>  a  and  v  "   0,  gives 

2    2 
P  (a  -  r  ) 

v  ■ 

4^1 
The  velocity  distribution  across  the  capillary  is  therefore  parabolic, 

The  volume  of  liquid  flowing  in  unit-tine  between  radii  r  and  r  +  dr  is 

3  -1 
2flrvdr  and  therefore  the  overall  flow  (  Q  en  sec  )  through  the  capillary 


is  given  by 


■/' 


2  Jfrvdr 


n*     (   (      l        l\A 
I  r(  a  -  r  ;dr 

»;i  J 


2ffP    1,2       2 
4^1 

77  Pa4 


87I 
This  is  known  as  Poiseuille  's  equation. 

Sources  of  Errors : 

1.  Streamline  and  turbulent  flow:  The  derivation  of  Poiseui lie's  equation 
assurr.es  that  flow  in  capillary  is  every  where  streamline.   It  is  found 
experiment ly  that  deviation  from  this  law  occur  at  high  rates  of  flow, 
this  was  shown  by  Osborne  Reynolds  (  1883  )  to  be  due  to  a  change  from 


streamline  to  turbulent  flow. 

2.  Kinetic  Energy  Correction:  Poiseuille's  Law  applies  strictly  only  to 
that  portion  of  the  tube  in  which  the  velocity  has  become  constant. 
In  most  types  of  viscometers,  part  of  the  applied  pressure  is  used  to 
give  kinetic  energy  to  the  liquid,  so  that  a  correction  must  be 
applied  in  order  to  obtain  the  pressure  used  in  overcoming  the  viscous 
resistance. 

3.  Viscous  End  Effects:   In  addition  to  the  work  done  in  overcoming  the 
viscous  resistance  in  the  capillary  itself,  a  small  amount  of  energy 
is  expended  in  overcoming  the  viscous  forces  between  the  converging 
diverging  streamlines  at  the  entrance  and  exit  respectively  of  the 
capillary. 

Couette  was  the  first  to  suggest  a  correction  to  account  for 
entrance   loss  for  a  Newtonian  fluid  in  a  capillary,  and  it  is 
generally  known  by  his  name. 

4.  Hydrostatic  Head  Correction:   In  viscometers  involving  applied  pressure 
the  effective  pressure  drop  across  the  capillary  is  the  sum  of  the 
externally  applied  pressure  and  due  to  the  head  of  liquid  in  the 
viscometer,  the  latter  of  course  normally  decreasing  during  the 
experiment  and  for  which  decrease  correction  must  be  made. 

5.  Drainage  Errors:  When  a  liquid  is  drained  from  the  bulb,  some  portion 
of  the  liquid  adheres  to  its  walls.  Several  types  of  viscometers 
include  a  second  bulb  of  similar  shape  above  the  main  bulb  to  make 
drainage  conditions  the  same  when  the  level  of  liquid  passes  the  upper  & 
lower  marks , 

6.  Surface-Tons  ion  Correction:   In  viscometers  in  which  there  is  free 
discharge  of  the  liquid  from  the  capillary  into  air,  the  formation  of 


LO 


droplets  is  associated  with  a  snail  pressure  drop,  and  tho  rate  of 
flow  for  the  sane  total  pressure  is  slightly  decreased  by  an  amount 
dependent  on  the  surface  tension  of  the  liquid.   This  effect  is 
eliminated  in  most  viscometers  by  using  submerged  discharge  of  the 
liquid  from  the  capillary. 
7.  Wall  Effect:  Analysis  of  capillary- viscometer  data  is  based  on  the 

assumption  that  the  velocity  is  zero  at  the  wall  of  the  capillary.   If 
a  fluid  has  a  finite  velocity  at  the  wall,  this  change  is  boundary 
condition  should  be  taken  into  account  in  deriving  viscometry  equation. 

b.   Glass  Capillary  Viscometers: 

Most  glass  capillary  viscometers  are  operated  by  the  force  of  gravity 
only  because  of  the  small  driving  force,  this  class  of  devices  is  useful 
for  low-viscosity  liquids  ranging  from  0.9  to  16,000  centistokes.   The 
viscometers  can  often  be  operated  by  application  of  external  pressure  in 
addition  to  the  hydrostatic  head.   In  this  way,  the  range  of  the 
viscometers  can  be  increased  considerably,  perhaps  as  much  as  40  times. 

The  principal  of  these  instruments  is  derived  from  the  viscometer 

3 

origionally  used  by  Ostwald  .   Basically,  the  viscometer  consists  of 

reservoir  bulbs  and  a  capillary  in  a  u-tube  arrangement,  as  shown  in  the 
Fig.  5.  The  efflux  time  of  a  fixed  volume  of  liquid  under  an  exactly 
reproducible  mean  hydrostatic  head  is  measured. 

This  origional  Ostwald  Viscometer  has  been  modified  in  many  ways  to 
minimize  certain  undesirable  effects  in  viscosity  measurements,  to  increase 
the  range  of -viscosity ,  to  meet  the  specific  requirements  of  certain  test 
liquids.   For  instance,  Caxnon  and  Fenske  modified  the  Ostwald  Viscometer 
so  that  the  upper  and  lower  bulbs  be  on  the  same  vertical  axis  in  order  to 
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Fig.   5.   OstwalcJ  Viscometer- 
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reduce  the  error  in  the  mean  head  caused  by  deviation  of  the  viscometer 
from  the  vertical  position.   For  opaque  fluids,  the  viscometer  was  modified 
to  have  reverse  flow  so  that  movement  of  the  meniscus  can  be  clearly 
observed . 

The  hydrostatic  head  available  for  liquid,  flow  in  a  glass  capillary 
viscometer  varies  with  time.   The  degree  of  hydrostatic  head  variation 
during  the  time  of  efflux  in  a  measurement  depends  upon  the  design  of  the 
viscometer.   For  a  Newtonian  liquid,  the  viscosity  of  which  is  constant, 
variations  in  the  pressure  head  have  no  effect  on  the  measurement, 
irrespective  of  the  fluid,  as  long  as  a  constant  volume  of  the  test  liquid 
is  used.   For  a  non-Newtonian  liquid,  it  differs  when  going  from  one  liquid 
to  another. 

Glass  Capillary  Viscometers  have  been  used  for  determining  viscosities 
of  Newtonian  liquids  because  of  their  excellent  accuracy,  relative 
cheapness,  and  simple  operation. 

Measurement  By  Gravity: 

Operation:  Although  specific  details  of  operation  vary  considerably  for 

the  viscous  types  of  glass  capillary  viscometers,  the  general  operating 

procedure  is  the  same.   Fig.  5  shows  a  schematic  diagram  of  an  Ostwald 

3 

Viscometer.   Typical  steps  in  making  a  measurement  are  (1)   a  sample  of 

liquid  to  be  tested  is  charged  through  tube  1  to  bulb  C;  (2)   the  viscometer 
with  the  sample  is  immersed  in  a  constant-temperature  bath  to  attain  the 
desired  temperature;  (3)   the  liquid  is  raised  in  tube  2  by  applying 
suction  until  the  level  is  above  mark  A,  and  by  removing  the  suction  the 
liquid  starts  to  flow  through  the  capillary  E;  (4)  the  efflux  time  of  the 
volume  of  bulb  D  between  the  A  and  3  marks  is  measured. 
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In  order  to  obtain  an  accurate  and  reproducible  result,  the 
hydrostatic  head  should  be  exactly  the  same  for  all  measurements.  This 

necessitates,  for  most  viscometers,  always  charging  the  instrument  with  a 

3 
given  amount  of  fluid.  The  kinematic  viscosity  of  the  liquid  is  calculated 

from  the  measured  efflux  time,  using  the  equation  given  below. 

V°  V?   «  k0  -  K/em  (.6   is  efflux  time  ) 

vhere  the  value  of  u-2   for  capillaries  with  trumpet-shaped  ends; 
m=l ,  for  square  ends,  and  in  general  it  is  assumed  that  m=l  in  relative 
viscosity  measurements  with  capillary  viscometers, 

K  is  assumed  to  vary  in  proportion  to  the  square  root  of  the  Reynolds 
number . 

Calibration: 

There  are  two  procedures  available 'for  calibration. 
(1).   By  Means  of  Standard  Fluids:  The  viscometer  is  calibrated  by 

measuring  the  efflux  time  of  a  standard  fluid.   In  order  to  determine 

the  two  coefficients  in  the  equation 

v-  1/p  -  k6   -  K/0m 

assuming  that  m=2 ,  two  measurements  are  made  with  oils  differing  by 
approximately  fivefold  in  viscosity.  The  coefficient  k  &  K  are 
calculated  for  m=2  by  solving  simultaneous  equations  obtained  by 
subsituting  viscosities  and  efflux  times  of  the  reference  materials 
into  the  above  equation. 
(2).  "Step-up"  Procedure:  '   In  cases  in  which  the  kinetic  energy 

correction  is  negligible  for  a  set  of  viscometers  (e.g.,  a  large 
length-to-diameter  ratio  for  the  capillary),  the  calibration  is  done 
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successively,  starting  with  water  as  the-  reference  fluid  in  the 
viscometer  with  the  smallest-diameter  capillary.  First,  the  efflux 
time  of  freshly  distilled  water  is  determined  with  smallest- 
capillary. 

Assuming  a  negligible  kinetic-energy  correction,  the  constant 
k  is  determined  in  the  above  equation  from  the  measured  time  with 
water.   Using  this  value  of  k,  the  viscosity  of  a  more  viscous  liquid 
is  determined.  This  viscosity  value  is  then  used  in  calibrating  a 
second  viscometer  with  larger  capillary  diameter.  After  the  second 
has  been  calibrated,  another  kinematic  value  can  be  determined  for  a 
previously  unstandardized  oil  of  still  higher  viscosity.  This  value 
is  than  used  to  calibrate  the  third  viscometer  and  so  on.  Steps 
between  successive  calibration  constants  or  viscosities  increase  by 
a  factor  of  three  or  less  until  the  desired  viscosity  range  is 
covered . 

3 

Orifice  Viscometers: 

In  the  course  of  development,  of  certain  process  industries,  the 
problems  of  measuring  and  controlling  flow  properties  of  various  classes 
of  fluids  have  been  encountered.  This  has  led  to  the  development  of 
instrument  for  resolving  these  problems,  and  a  number  of  them  have  come 
into  general  use.   For  instance,  the  oil  industry  adopted  the  method  of 
measuring  viscosity  with  short-tube  or  orifice  viscometer.  These 
instruments  are  now  common  in  industry:  e.g.,  the  Saybolt-universal  and 
Furol  in  the  United  States ; Redwood  No.  1  and  Ko.  2  in  Great  Britain, 
Englsr  in  Germany,  and  Barbey  in  France.   In  the  paint  industry  cup-type 
viscometers  have  been  developed  and  are  widely  employed  at  the  present  time. 
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The  origional  design  concepts  of  those  viscometers  were  derived  from 
the  Hagen-Poiseuille  Law  which  states  that  the  efflux  tine  of  a  fixed 
volume  through  a  capillary  is  proportional  to  the  viscosity  of  the  fluid. 
Unfortunately,  the  viscometers  that  were  developed  consisted  of  short 
capillaries  or  orifices.   Flow  in  these  viscometers  does  not  obey  the 
Hagen-Poiseuille  Law  and  efflux  times  are  not  in  any  simple  relation  to 
the  viscosities.   An  added  factor  in  the  development  of  such  viscometers 
was  the  practical  requirement  that  the  method  should  be  simple,  quick  and 
reliable. 

In  these  viscometers  the  orifice  length  does  not  exceed  10  ti-mes  the 
diameter  of  the  orifice.  The  hydrostatic  head,  i.e.,  the  driving  force 
causing  fluid  to  flow,  is  consumed  to  an  appreciable  extent  at  the  orifice 
enterence.   The  friction  loss  due  to  this  effect  is  a  function  of  the 
cross-sectional  area  ratio  of  cup-to-orifice,  velocity  of  the  fluid,  shape 
of  the  orif ice-enterence  port  of  the  viscometer,  and  properties  of  the 
fluid.   Furthermore,  the  situation  becomes  more  complicated  because  of  the 
varying  hydrostatic  head  during  an  experimental  run.   For  this  reason,  the 
efflux  time  readings  for  these  viscometers  are  not  in  the  same  relation  to 
the  driving  force  as  with  the  usual  capillary  viscometers,  and  conversion 
formulas  or  tables  must  be  consulted  in  comparing  results. 

This  type  of  viscometer,  however,  is  widely  used  because  of  simplicity 
and  inexpensive  operation.  These  viscometers  are  mostly  used  for 
Newtonian  fluids. 

Principle  of  Operation: 

<  Viscometers  of  this  type  consists  of  a  reservoir  and  orifice  (  about 
10  diameters  long  at  most  )  with  or  without  a  temperature-control  jacket 
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and  a  receiving  flask.  This  method  of  operation  is  nearly  the  same  for 

all. 

(1).  The  liquid  under  test  is  poured  into  a  cup  sorrounded  by  a  water  or 

oil  bath  providing  temperature  control. 
(2).  The  level  of  the  liquid  in  the  cup  is  adjusted  to  a  definite  height. 
(3).  When  the  desired  temperature  is  attained,  the  valve  at  the  base  of 

the  cup  is  opened. 
(4).  The  time  required  for  a  specified  volume  of  liquid  to  discharge 

through  orifice  into  a  measuring  vessel  placed  below  is  measured. 
The  measured  efflux  time,  generally  in  seconds,  is  a  purely  arbitrary 
expression  of  the  viscosity.   In  order  to  convert  these  figures  into 
absolute  units,  an  empirical  formula  for  each  instrument  is  derived.  The 
equation  is  usually  in  the  form  of 

/ 

This  form  is  the  same  as  that  for  glass  capillary  viscometers  with 
the  difference  that  an  appreciable  end  effect  is  reflected  in  k.  So  it 
requires  calibration  for  various  ranges  of  efflux  time.   For  approximate 
values  viscosity  conversion  charts  are  available  for  reference  for  orifice 
viscometers. 

Rotational  Viscometers: 

Sir  Isaac  Kewton  first  observed  the  relation  between  velocity  and 
resistance  to  flow  for  a  rotating  solid  cylinder  in  a  uniform  and  infinite 
fluid.  He  showed  the  exponential  decrease  in  rate  of  shear  for  the  rotating 
concentric  cylindrical   laminar  shells  of  fluid  as  one  moves  away  from  the 
cylinder  and  into  the  fluid.   In  fact  a  complete  solution  to  Margule's 
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equation  was  not  arrived  at  and  it  was  Store's  who  discovered  it  later. 

Nearly  two  hundred  years  elapsed  before  the  first  practical  rotational 
viscometer  was  devised  by  Couette  in  1890.  Couette's  concentric-cylinder 
viscometer  consisted  of  a  rotating  cup  and  an  inner  cylinder  which  was 
supported  by  a  torsion  wire  and  rested  in  a  point  bearing  in  the  bottom  of 
the  cup.   In  this  unit  the  principle  of  guard  rings  was  utilized  to 
eliminate  end  effects  as  shown  in  the  Fig.  6.  The  guard  rings  F  and  F» 
are  fixed  and  do  not  move,  where  as  the  inner  cylinder,  represented  by  A, 
turns  in  proportion  to  the  viscosity  of  the  fluid  in  the  gap  until  the 
opposing  forces  of  the  torsion  wire  come  in  balance  and  there  is  no  net 
movement  of  the  inner  cylinder.  When  the  system  is  thus  balanced  the  fluid 
in  the  gap  flows  by  the  stationary  inner  cylinder  and  the  guard  rings. 
There  is  no  flow  over  the  top  and  bottom  of  the  inner  cylinder  since  the 
fluid  inside  thy  guard  rings  is  trapped  ,and  is  stationary.   Couette's 
design  enabled  him  to  calculate  the  apparent  viscosities  of  non-Newtonian 
samples  with  only  a  small  error  because  of  the  very  small  ratio  of  gap  to 
inner  radius. 

In  1913  Hatscheck  described  a  modified  version  of  the  Couette 
viscometer  which  improved  on  the  design  of  the  guard  rings  and  minimized 
effects  with  very  low  viscosity  fluids.  Until  1940  there  was  not  much 
literature  on  rotational  viscometery.  Since  then  a  good  many  designs  have 
appeared,  some  of  which  are  available  commercially. 

Theory : 

A  rotating  body,  immersed  in  a  liquid,  experiences  a  viscous  drag  or 
retarding  force.  The  amount  of  viscous  drag  is  a  function  of  the  speed  of 
rotation  of  the  body. 


18 


Torsion  wire 


Guard  rir.2 
support 


» J- 


'     Perforation  for^ 
ring  support 


L 


Fig.    6.    Couette's   Viscometer" 
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Rotational  viscometers  arc  usually  more  elaborate  than  the  capillary 
type.  When  used  with  Newtonian  fluids  they  are  often  less  accurate,  but 
they  have  a  number  of  advantages  including:  (1)  measurements  are  made 
under  steady  state  conditions;  (2)  during  a  given  measurement  variation  in 
shearing  stress  throughout  the  material  may  be  made  quite  small;  (3)  the 
same  sample  may  be  used  for  measurement  at  different  shear  rates  or  for 
continuous  measurements  on  materials  whose  properties  change  with  time; 
(4)  measurements  made  on  non-Newtonian  materials  can  be  interpreted  in 
absolute  units.   For  these  reasons  rotational  viscometers  have  become  the 
single  most  widely  used  class  of  instruments  for  Theological  determinations. 

Derivation  of  Basic  Equation: 

Although  rotational  viscometry  includes  spheres,  disks,  cones  and 
odd-shaped  rotors  and  cups,  the  most  common  type  is  the  coaxial-cylinder 
viscometer.   In  its  simplest  form  as  shown  in  the  Fig.  7,  the  coaxial 
cylinder  viscometer  comprises  an  inner  cylinder  (  radius  R,cm,  height  h  cm  ) 
and  an  outer  cylinder  (  radius  R2cm  ).  The  outer  cylinder,  which  acts  as 
the  container  for  the  liquid  under  test,  is  rotated  at  a  constant  speed 
(/i-rad  sec  ),  and  the  resultant  torque  (  T  dyn  cm  )  is  measured  by  the 
angular  deflection  of  the  inner  cylinder  which  is  suspended  by  a  fine  wire. 

The  following  assumptions  are  made  in  order  to  arrive  at  the 

4  9 
fundamental  equation:'  (1)  the  liquid  is  incompressible;  (2)  the  motion  of 

the  liquid  is  laminar;  (3)  the  streamline  of  flow  are  circles  on  the 

horizontal  plane  perpendicular  to  the  axis  of  rotation  (  i.e.,  the  velocity 

is  a  function  only  of  radius jradial  and  axial  flows  are  assumed  to  be  equal 

to  zero  );  (4)  the  motion  is  steady-all  time  derivatives  in  the  equations 

of  continuity  and  motion  are  zero;  (5)  there  is  no  relative  motion  between 
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Fig.    7.    Co-axial  Cylinder  Viscometer 
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the  surface  of  the  cylinders  and  the  fluid  in  immediate  contact  with 
cylinders,  i.e.,  no  slippage;  (6)  the  motion  is  two-dimensional 
(  neglecting  the  edge  and  end  effects  ). 

Let  V   be  the  angular  velocity  of  the  liquid  at  a  distance  r  from  the 
axis  of  rotation,  and  let  the  shearing  stress  and  rate  of  shear  at  this 
radius  be  Fr  and  Dr  respectively.  Considering  the  part  of  the  liquid 
between  the  inner  cylinder  and  an  arbitrary  radius  r.   The  torque  exerted 
on  this  liquid  by  the  outer  layers  of  liquid  is 

shearing  stress  x  surface  area  x  radius  ■  Frx  2  jj  rh  x  r 

=  2rfr2h?r 

The  torque  actixig  on  the  inner  surface  is  the  measured  torque  T. 
Since  the  motion  is  assumed  steady,  these  torques  must  be  equal 
therefore 


T  =  2  ffr2hFr 


■ 


F. 


2  tfr2h 

Thus  the  shearing  stress  is  inversely  proportional  to  the  square  of 
the  distance  from  the  axis  of  rotation.  Now  by  making  the  annular  gap 
small  as  compared  to  the  cylinder  radii,  it  is  observed  that  the  shearing 
stress  may  be  made  almost  constant  throughout  the  material  being  sheared. 
This  is  in  contrast  to  the  capillary  viscometer  for  which  the  shearing 
stress  varies  from  zero  at  the  axis  to  a  maximum  at  the  wall  of  the 
capillary. 

To  find  an  expression  for  the  shear  rate,    differentiate  the  linear 
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velocity  with  respect  to  r. 
So 


v  ■  rw 


dr      dr 
The  term  r  (-—)  represents  the  shear  rate  since  the  other  term,  w,  is 
the  radial  velocity  gradient  of  a  rigid  rotating  object. 
Thus 

n      dw 
D_.  =  r  — 

r     dr 


For  a  Newtonian  liquid 


F  «  D  V  mid  therefore 

r    r  fc 


dw  ^ 


2/Tr2h   '  dr 


dw  m  __T 

dr    2//»j>hr3 

Integrating  this  using  the  no-slip  boundary  condition  that 


w  ■  0   when   r  ■=  R* 
and  w  -  Xi-    when   r  =  R„        find  that 

T  -  agRlWhYA,  c  ^  /j. 
R2  -  RX 

Where  C  is  an  instrument  Constant. 

The  viscosity  is  therefore  directly  proportional  to  the  ratio  of 
torque  to  angular  velocity.   It  may  be  noted  that  the  above  equation  is 
equally  applicable  when  the  inner  cylinder  is  caused  to  rotate  at  angular 
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velocity  jO- and  the  outer  cylinder  is  at  rest. 

End  Effects: 

All  of  the  preceeding  mathematical  derivations  assumed  characteristics 
of  an  infinitely  long  cylinder  in  that  no  account  was  taken  of  the  drag  on 
the  ends  of  the  inner  cylinder.  Where  the  term  h  was  used  for  the  height, 
the  expression  h  +  hQ  should  have  been  used  where  hQ  is  the  additional 
height  to  be  added  due  to  ends.   This  effect  is  taken  into  account  in  the 
single  instrument  constant  C  which  is  determined  by  calibration  with  one 
or  more  liquids  of  known  viscosity. 

Temperature  Effects : 

To  a  greater  or  lesser  extent,  the  viscosity  coefficient  of  all  fluids 
are  dependent  on  temperature.  So  thermostatic  control  of  the  material 
during  measurement  is  highly  desirable.  While  this  presents  no  great 
difficulty  with  const ant -tor que  type  viscometers  in  which  the  cup  is 
stationary,  it  is  less  easy  with  rotating  cup  instrument. 

Friction: 

Viscosity  measurements  involving  small  torque  values  are  often 
difficult  because  of  appreciable  friction  arising  in  the  bearings.  While 
jewelled  bearings  have  been  used  successfully,  one  of  the  most  satisfactory 
method  is  that  described  by  Oldroyd ,  Strawhridge  and  Toms  (  1951  )  in 
which  the  main  spindle  carrying  the  bob  is  positioned  by  a  number  of  very 
small  jets  of  air  directed  radially  on  to  it. 

Wall  Effects: 

"While  deriving  the  general  equation  it  was  assumed  that  the  bulk 
viscous  properties  of  the  material  apply  right  up  to  the  surface  and  that 
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i 
no  slippage  occurs  at  any  solid  boundary.  Wall  effects  nay  be  suspected  if 
the  (F,D)  curves  determined  with  different  annular  gaps  differ  appreciably. 

Time-Dependent  Effects  (  Thixotropy  ): 

There  are  reversible  and  irreversible  effects  varying  with  time.  Of 
these  reversible  work  softening  or  thixotropy  is  most  frequently  discussed 
in  the  literature. 

The  rotational  viscometer  is  a  particularly  suitable  instrument  for 
the  study  of  thixotropic  materials  since  it  permits  measurements  to  be 
made  continuously  on  the  same  sample  of  material  under  varying  shear 
conditions.  There  are  at  least  two  approaches  to  the  semiquantitative 
determination  of  thlxotropic  changes.  These  are  (1)  the  measurement  of 
the  hysteresis  between  "up"  and  "down  curves,"  according  to  the  method  of 
Greenland  (2)  determination  of  the  decay  of  shear  stress  as  function  of 
time  at  one  or  more  constant  shearing  rates  (  Green  and  Weltman  ). 

Conclusion: 

Rotational  viscometers  are  versatile  laboratory  instruments  for 
measuring  flow  properties  of  many  fluids.  The  simplicity  of  design  and 
ease  of  manipulation  of  the  instrument  has  great  appeal  to  most 
reheologists.  When  treating  data  obtained  in  the  rotational  viscometer, 
it  is  important  to  keep  in  mind  the  effect  of  deviations  from  Newtonian 
behaviour.   The  rate  of  shear  is  not  uniform  across  the  gap  and  the 
Newtonian  expression  cannot  be  applied  directly.  Shear  stress  and  the 
rate  of  shear  should  be  expressed  at  the  same  point  in  the  viscometer  and 
the  corrections  discussed  earlier  should  be  taken  into  account,  the  most 
important  of  which  is  the  relation  of  gap  size  to  radius  of  inner  cylinder. 
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Calibration  of  the  instrument  should  be  done  carefully  and  other  factors 
liice  temperature  turbulence  and  end  effects  should  also  be  considered. 

Falling  Ball  Viscometers: 

The  flow  behaviour  of  a  fluid  in  a  rheological  instrument  is  obtained 
in  principle  by  solving  the  equation  of  notion- for  the  given  geometry  -with 
appropriate  boundary  conditions.  When  a  solid  body  is  allowed  to  fall 
under  gravity  through  a  viscous  medium  a  period  of  initial  acceleration  is 
followed  by  motion  at  uniform  terminal  velocity  for  which  the  gravitational 
force  is  balanced  by  the  viscous  resistance.  Measurement  of  this  terminal 
velocity  affords  a  means  of  determining  the  viscosity  of  the  medium. 
Bodies  of  various  shapes  may  be  used;  the  simplest,  both  from  theoretical 
and  experimental  points  of  view,  is  the  sphere.  This  method  however  is 

limited  for  use  to  Newtonian  fluids. 

i 
Derivation  of  the  Basic  Equation: 

The  principle  of  the  falling-sphere  viscometer  is  the  well-known 

Stolce's  Law  which  relates  the  viscosity  of  a  Newtonian  fluid  to  the 

velocity  of  the  falling  sphere.  According  to  Stol:e»s  the  viscous  resistance 

to  the  motion  of  a  sphere  moving  with  velocity  v  is 

6jyr»jv 

Where  r  is  the  radius  of  the  sphere.  The  driving  force,  due  to  the 
difference  in  density  between  sphere  and  the  fluid,  is 


|rfr3(cr-  ?)< 


Where 

cr  c  density  of  sphere 


i 

?  -   density  of  fluid 
g  «=  acceleration  due  to  gravity 
Equating  the  two  expressions,  the  viscosity 

„  2gr2(v-  P) 
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Which  can  be  known  by  measuring  velocity  (v) ,  r.crand  f  being  fixed. 
In  the  derivation  of  the  above  equation  it  was  assumed  that  the 
sphere  moves  with  a  very  slow  velocity  in  the  fluid  of  infinite  extent. 
Since  this  ideal  situation  is  not  realised  in  actual  viscometers,  the 
following  corrections  should  be  taken  into  consideration  in  calculating 
the  absolute  viscosity. 

(1).  The  motion  of  the  sphere  relative  to  the  liquid  is  slow  (  Finite 
Reynolds  number  ) : 

It  was  assumed  in  Stoke*s  Law  that  the  Reynolds  number  given  by 

KR  «=  2  rvP/q  is  smaller  than  unity. 

3 

To  overcome  this  restriction  to  a  rate  of  flow  Cseen  and 

Goldstein3  developed  a  correction  factor   in  terms   of   infinite  series, 
as  gi Veil  below: 


?abs=    VC  l  +  (3/l6)l^e-    (19/l280)KRe2   +   (7l/20480)NRe" 


(2). 


Wall  Effects: 

It  is  obvious  that  in  any  experimental  determination  of 
viscosity  the  liquid  cannot  be  of  infinite  extent.   For  a  sphere 
falling  in  a  cylinder  of  finite  dimensions,  the  force-resisting  motion 
of  the  sphere  is  found  to  be  increased  by  the  drag  exerted  by  the 


27 


«3 

cylinder  wall.   ladenburgJshovred  that  the  Stoke's  viscosity  should  be 
corrected  as  follows: 

.?abs  "  */f  l+2.1r/R  J 

Where  R  is  the  radius  of  the  container. 

(3).  Homogenity: 

It  was  assumed  that  the  liquid  should  be  homogeneous.   But  many 
liquids  change  their  viscosity  with  age  and  others  show  an  apparent 
viscosity  which  varies  with  rate  of  shear. 

(A).  End  Effects: 

According  to  Ladenburg,  the  resisting  force  is  increased  by  the 
factor  as  given: 

7abs=  V/  l+3.3r/h  J 

Where  h  is  the  height  of  the  sphere  from  the  bottom  of  the  cylinder. 

Falling-ball  measurements  are  generally  employed  for  reasonably 
viscous  materials  because  of  the  difficulty  of  handling  the  very  small 
balls  or  of  correctly  measuring  the  small  differences  in  density  between 
ball  and  liquid  needed  to  get  a  suitably  slow  rate  of  flow  in  a  highly 
fluid  medium.  This  technique  is  also  confined  to  cases  in  which  disturbance 
of  the  fluid  prior  to  the  viscosity  measurement  is  to  be  avoided  that  rule 
out  the  use  of  other  types  of  viscometers. 
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CHAPTER  II 
EXPERIMENTAL  PROGRAM 

General: 

In  fact  before  starting  the  actual  experiment  with  the  Indentor 
Viscometer,  the  viscosities  of  the  three  Test  Fluids  were  determined  with 
the  Brookfield  Synchro-Lectric  Viscometer.  2  A  brief  description  of  the  same 
is  given  below: 

1.  Brookfield  Synchro-Lectric  Viscometer  Method: 

This  method  is  comprised  of  the  following  steps: 

(a).  Test  Samples: 

These  sample  fluids  were  obtained  from  Brookfield  Engineering 
Laboratories  Incorporated,  Stoughton,  Masschusettes  ,  U.S.A. 


Table 

No 

:   1 

Test  Fluids 

• 

Sr  No 

Sample 

Lot  No. 

Standard 
Viscosity 

Temperature 

1 

R2 

120567 

29,100  cps 

77°F 

2 

Hl 

080767 

59,400  cps 

77°F 

3 

H2 

121367 

95,000  cps 

77°F 

(b).   Experimental  Set  Up: 

This  consisted  of  the  following: 

•1.  Brookfield  Viscometer  (  RV  type  ) 

2.   Test  Fluids 
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3.  Three  600  c.c.  beakers 

4.  A   constant  temperature  bath  with  thermometer 

5.  Brookfield  Factor  Finder 

The  various  working  components  of  the  Brookfield  Synchro-Lectric 

12 

viscometer  are  as  shown  in  the  Fig.  8. 

(c).  Experimental  Procedure: 

1.  The  spindle  RV-7  (  which  only  covered  the  range  of  viscosity  to  be 
measured  )  was  attached  to  the  lower  shaft  while  the  shaft  was  held 
firmly  and  avoided  the  side  thrusts  in  order  to  protect  the  alignment. 

2.  The  spindle  was  then  lowered  into  the  test-fluid  avoiding  air 
bubbles  until  the  fluid  level  was  at  the  groove  cut  in  the  spindle 
shaft. 

3.  The  viscometer  was  then  levelled  with  the  help  of  the  air  bubble 
level  and  set  the  speed  of  rotation. 

4.  Then  the  clutch  was  depressed  and  started  the  viscometer  motor. 
On  releasing  the  clutch  the  dial  rotated  till  the  pointer 
stabilized  at  a  fixed  position  on  the  dial.  Different  readings 
were  taken  for  different  speeds  of  rotation  ranging  from  10-100r. 
p.m.  depending  upon  the  test  fluid  used.  At  high  speed  the  pointer 
was  held  on  to  the  dial  when  in  view  by  depressing  the  clutch. 
Then  the  motor  was  switched  off  and  the  reading  noted  when  the 
rotation  stopped.  The  readings  were  chocked  by  repeating  the 
procedure. 

5.  The  viscosity  of  the  test  material  was  then  found  by  multiplying 
the  dial  reading  and  the  mult ipicat ion  factor  found  by  consulting 
data  provided  by  the  Brookfield  Company. 
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(d).  Data  and  Calculations: 


Table-2 
Observations  for  Fluid  Sample-1  Using  Spindle  RV-7 


Multipication  Dial     Calculated  Standard   %age 
Sr  Ko  R.PM     .,  Remarks 

i-actor      Reading   Viscosity  Viscosity  Deviation 

1  10      4000        7.0    28,000  cps  29,100  cps  3.80 

Temp. 

2  20  2000  14.1  28,000  cps     29,100  cps     3.44  of 

bath 

3  50  800  35.4         28,320  cps     29,100  cps     2.68_  77°F 

4  100  400  70.9         28,360  cps     29,100  cps     2.55 


Table-3 
Observations   for  Fluid  Sample-2  Using  Spindle  RV-7 


e     „  ^,     Multipication     Dial  Calculated     Standard       %age 

Sr  iNo     R.PM  Remarks 

Factor  Reading       Viscosity     Viscosity     Deviation 

1  10  4000  14.3         57,200  cps     59,400  cps     3.7  Bath 

Temp . 

2  20  2000  28.7         57,400  cps     59,400  cps     3.37  77°F 

3  50  800  71.8         57,440  cps     59,400  cps     3.3 

out  of 

4  100  400  "  "  "  "  Range 
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Table-* 

Observations   for  Fluid  Saraple-3  Using  Spindle  RV-7 


Sr  No 

R.PM 

Multipicat 
Factor 

ion 

Dial 

Reading 

Calculated 
Viscosity 

Standard 
Viscosity 

95,000  cps 

%age 
Deviation 

Remarks 

1 

10 

4000 

22.8 

91,200  cps 

4.00 

rBath 

2 

20 

2000 

45.7 

91  ,400  cps 

95,000  cps 

3.79 

Temp. 
77°F 

3 

50 

800 

ii 

ti 

tt 

ti 

Out 

4 

100 

400 

ii 

tt 

« 

•i 

of 

Range 

33 


2.   Indentor  Viscometer  Method: 

This  can  be  divided  into  following  steps: 
(a).  General  Arrangement: 

In  this  method  the  same  three  Test  Fluid  Samples  were  used.   Besides 
this  there  were  three  indentor  sizes,  three  different  weights  and  three 
depths  of  indentation  for  computation  of  velocities  at  these  different 
depths  as  shown  in  table-5. 

Table-5 
Test  Fluids  and  Indentor  Sizes 


Test  Fluid    Indentor    Weight  of  Depth  of 

No:     Viscosity     Radius      Indentor  Indentation 

. (£1 l*u &i_ m . 

1  29,100  cps          0.75    Inches      103.5gms  0.05 Inches 

t 

2  59,400  cps    1.00  Inches   138.5gms  O.OSInches 

3  95,000  cps    1.25  Inches   173.5gms  O.lOInches 


Each  experimental  sequence  constituted  a  combination  of  one  indentor 
size,  one  weight  of  indentor,  one  depth  of  indentation  and  the  testing 
fluid  at  a  time  resulting  into  81  total  number  of  experiments  to  be 
performed.  The  various  combinations  of  these  sequences  are  shown  in  table-6, 

If  the  experimental  sequence  has  a  combination  as 

F  R  V?  D 
12  3  2 

*  it  means   that  the  fluid   used    is  no:l    (  29,100  cps  viscosity   ),   Indentor 

no:2    (   1.0"Radius    ),  weight   of   indentor    is   no:3    (   173.5   gms    )   and   depth  of 
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indentation  is  no:2  (  0.08"). 

The  experiments  were  performed  in  a  purely  random  manner.  This  vas 

achieved  by  consulting  the  "Random  Digits"  tables  from  the  book  The 

13 
Compleat  Strategy.? t  by  J.D.  Williams.   This  random  order  is  as  shown  in 

table-7. 


Table-7 

Random  Experimental  Order 


Exp  Seq  ExpSeq  Exp  Seq  Exp  Seq  Exp  Seq  Exp  Seq  Exp  Seq  Exp  Seq  Exp  Seq 
No  No  No  No  No  No  No  No  No  No  No  No  No  No  No  No  No  No 


1  11  10  44  19  56  23  67  37  45  46  09  55  14  64  06  73  41 

2  16  11  36  20  31  29  75  38  27  47  15  56  38  65  35  74  70 

3  43  12  79  21  28  30  60  39  54  48  24  57  19  66  25  75  07 

4  63  13  22  22  72  31  71  40  61  49  66  58  23  67  58  76  30 

5  18  14  62  23  64  32  32  41  57  50  34  59  03  68  53  77  02 

6  21  15  73  24  20  33  55  42  04  51  42  60  05  69  76  78  69 

7  59  16  40  25  37  34  52  43  78  52  29  61  46  70  74  79  80 

8  17  17  47  26  51  35  65  44  13  53  50  62  08  71  33  80  77 

9  10  18  49  27  4S  36  26  45  68  54  12  63  01  72  39  81  81 
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b.  Experimental  Set  Up: 

The  appratus  used  in  this  method  is  shown  in  the  Fig.  9  and  is  listed 
as  under: 

1.  Indentor  Viscometer  Unit:  This  consisted  of  the  following  (Fig.  10) 

a.  Supporting  Rack  with  levelling  base. 

b.  An  Electromagnet. 

c.  Displacement  transducer. 

d.  Guides  for  the  indentor. 

2.  Strip  chart  recorder 
Brush  Recorder  Mark  1 

3.  Four  1  I2   volts  batteries  for  supply  to  the  displacement  transducer. 
A.  Supply  unit  for  the  Electromagnet  NJS  corporation  Model  RB  36-2 

0-36  VDC  0-2  ,\rap. 

5.  Three  glass  containers  of  equa^.  size  90  mm  x  50  mm  and  their 
covers , 

6.  Spatula. 

7.  Kleenex  box. 

c.  Experimental   Procedure: 

1.  The  instrument  was  levelled  in  order  to  avoid  any  error  due  to 
misalignment. 

2.  According  to  the  experimental  sequence  the  required  fluid, 
indentor  and  weight  were  selected. 

3.  After  mounting  the  indentor  and  the  weight,  the  instrument  was 
calibrated  by  adjusting  the  travel  0.1"  from  the  one  end  of  the 
chart  on  the  recorder  to  the  other  end  covering  80.0  mm.   This 
calibration  was  checked  once  again  before  starting  the 

exper  imentat  i  on . 
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Fig.   10     Indentor  Unit 
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4.  The  indentor  was  then  adjusted  to  just  touch  the  fluid  surface. 

5.  The  chart  speed  was  set  on  200  nun/sec  and  the  indentor  released 
from  the  magnet  by  pressing  the  switch. 

6.  The  velocity  profile  was  observed  on  the  chart  and  the  recorder 
switched  off  on  completion  of  the  travel. 

In  this  manner  all  the  81  experiments  were  performed  under  similar 
conditions  of  temperature  and  calibrating  the  instrument  every  time  in 
order  to  avoid  any  variations.  Moreover  the  travel  of  the  indentor 
disk  was  kept  0.1"  throughout  but  measurements  were  made  according  to 
the  required  depths  from  the  chart, 
d.  Data  and  calculations  are  aontained  in  Appendix  B. 
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CHAPTER  III 


ANALYSIS 

Method  Of  Least  Squares: 

The  analysis  of  the  data  is  "based  upon  the  derivation  of  the  Indentation 

Ik 

of  a  Newtonian  Fluid  by  a  Right  circular  cylinder. 

F, 


J 


v>  *  "8  *j  "} 


where  j  is  the  experiment  number 


V   =  Coefficient  of  viscosity  of  the  fluid 
F  =  Force  on  the  indentor 
a  =  Radius  of  the  indentor 
v  =  Velocity  of  indentation 


It  is  assumed  that  viscosity  may  be  dependent  upon  the  fluid  used, 

2       2       2 

v,   a,  F,  v   ,   a  ,   F    ,   av,   aF,   vF  and  some  experimental  error  e. 

Then 

t 

^=CILlj+C2L2J+SSj+C^j+SFj+C6Vj+SVj2+C8aj2+SFj2+C10ajrj+CllajVj+C12V/j+eJ 


*UjHLlJ   L2J   L3J    aj   FJ   Vj   V/   a/  F/   Vj    Vj   Yj)    f     ll 


*H 


K  C12> 


This  represents  a  linear  model   and  in  the  matrix  form  it  becomes 


rCli 


l»j|-  (  A  }  ' 


'e,   -> 


>  +  < 


12j       Le8l 


>> 
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V   =  (  A  )  {C}  +  {e} 


where 


p  is  the  dependent  response  variables 
A  is  the  input  independent  variables 

If  errors  are  present  in  the  dependent  variables  and  independent 

variables  then  the  problem  becomes 

{p  +  6}  ---  (  A  +  A)  {C}  +  {e} (1) 

in  which  {6}  is  a  (  8lxl  )  column  matrix  of  unknown  errors  in  y, 

(A)  is  a  (  81x12  )  rectangular  matrix  of  unknown  errors  in  (  A  ) . 
Multiplying  (l)  by  (  A  +  A),  yields 

(A  +  A)  {p+6}  =  (A  +  A)  (A  +  A)  {C}  +  (A  +  A)  {e}  (2) 

A 
A 

The  estimate  {C}  which  solves  the  equation  (2)  in  the  least  square 
sense  is  that  for  which 

(A  +  A)  {e}  =  0  , 

substituting  this  in   (2)   and  solving 

{C}  =    [(A  +   A)    (A  +  A)]"1   (A  +  A)    Ui+6} 

{C}   =   (A  A+AA+AA+AA)"1   (A  +  A)    {y+6} 
saving  the   first  order  terms 

{C}  =   (A  A+AA+AA )_1   (A  +  A)    {y+6}      (3) 

Now  let 

(A  A+AA+AA)-1  =   (M+e)-1       where  M  =  A  A,   e  =  AA+AA 

(AA+AA+AA)"1  =   (M   (I+M~1e))~1 
=   (l+M~1e)~:LM'':L 
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By  binomial  expension  and  saving  the  first  order  terns  this  becomes 

(AA+AA+AA)"1  =  (I-M_1e)  M*"1 
Therefore  (3)  becomes 

{C}  =  (I-(A  A)_1(AA+AA))  (A  A)_1(A+A)  {y+6} 

On  multiplication  and  leaving  the  terms  other  than  the  first  order, 
yields 

{C}  =   [(A  A)""2*  +(A  A)_1A  -   (A  A)_1(AA+AA)(A  A)-1Aj{y+6} 

{C>   =   (A  A)-1A  p+(A  A)_1A  u-(A  A)~1(AA+AA)  (A  A^Ay+U  A)-1A6 (k) 

The  first  term  in  equation  (k)   is  the  least  squares  estimate  of  the 
coefficients  obtained  when  errors  are  present  only  in  the  dependent  variables 
and  the  model;  the  second  and  third  terms  give  the  effects  of  errors  in  the 
independent  variables,  with  the  second  term  indicating  the  combined  effects 
of  ill-conditionedness  of  the  experiment  and  these  errors;  the  fourth  term 
the  effect  of  errors  in  the  dependent  variables  on  the  estimate  of  the 
coefficients . 

Eq'n  {h)   can  be  used  to  find  the  confidence  intervals  for  the 
estimates  of  the  coefficients  {C}. 
Therefore 

Sj  ■  V  Vj A  +  Ei ViAk 


*   2  2  2  2  2 

c  Jk       6k  Ik       A.k 

2 

The  values  of  o"   can  be  used  to  estimate  confidence  intervals  for  the 

cj 

parameters  C  as 

C  =  C  +  a: 
J    J    p 
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CHAPTER   IV 
RESULTS    AND  COiJCWSIONS 


Computed  Results: 

The  analysis   of  the  equation 


MS   -  C1Llj+C2l^j+C3L3j+C4aj+C5Fj+C6vj+C7vj2+C8aj^C9Fj2+C10ajFj+C11a.vj+C12vjFj 


has  revealed  the   following  results: 


Estimate  of  the 
coefficient 


Ct  -    .472 

C2  -    .353 


C3 


C5 
C6 
C7 

S 
C10 
Cll 
C12 


> 


a 

.255 

+ 

.295 

- 

.004 

+ 

.004 

- 

.0000006 

- 

.026 

+ 

. 000004 

+ 

.001 

— 

.001 

+ 

.000001 
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Table-8 
The  Standard  Deviations  for  the  Estimates 


of  the  Co-e: 

Eficients 

*i* 

ID"' 

8.748 

0.046  x 

9.166 

0.048  x 

-1 

10 

8.857 

0.045  x 

-1 
10 

5.825 

0.024  x 

-1 

10 

0.022  x 

-2 

10 

0.002  x 

-4 

10 

0.485  x 

-3 

10 

0.017  x 

-5 
10 

0.147  x 

-8 
10 

0.648  x 

-12 

10 

0.2401 

0.115  x 

-3 

10 

0.605  x 

-9 

10 

-11 
0.2404  x  10 

0.3324 

-4 
x  10 

-7 
0.2788  x  10 

0.7495 

,     -4 
x  10 

-7 
0.3188  x  10 

0.2538 

-8 
x  10 

-11 
0.1865  x  10 

Cl 

C2 
C3 
C4 
C5 
C6 

C7 

C 

8 

C9 
C10 
Cll 
C12 


For  further  details  of  results  appendix  B  may  be  referred. 
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Discussion  And  Conclusions: 

On  observing  the  computed  values  of  the  coefficients  C^ ,  Cg,  C-j ,  Cq  , 
C, 0,  C,,  and  C-2  which  can  be  neglected  being  very  snail  as  compared  to 
other  coefficients,  the  resulting  equation  is  written  as 

ui  "  ciLij .+  P2*2j  +  cihi  +  c4a  4  c8a2  * 

«  -.472  Ljj   -   .353  1^.   -   .255  Lj.  +   .295  a  -   .026  a2 
As  the  radii    of  the  indentors   lie  between 

1.905  ^  a  ^  3.175 

So 

2 

U    -  -.472L1j-.353L2.-.255L3,+.295(a-1.905)+.5619-.026(a-1.905)   -.02   (3.S10) 

(a-1. 905)-. 026(3. 629) 
-  -.472L1.-.353L2j-.255L3.  +  .295(a-1.905)-.026(a-1.905)  -.101  (a-1 .905)-!-. 4676 


-   (-. 472+. 4676)Lj.+(-. 353+0. 4676)L2-  +  (-. 255+0. 4676)L3j-i-0.194(a-l. 905)-. 026 

(a-1. 905) 

or 

2 

A        sec  «=  (-.0044)L1.+(0.1146)L2.+(0.2126)L3.+0.194(a-1.905)-.026(a-1.905) 

Converting  this  equation  of  jj.  with  units  as  Centipoises  by  multiplying 
vlth  a  factor  of  (  980.6  x  100  ),  it  becomes 

Centipoises  "  -WlLjj+ll  ,23312.+20, 8451^+19, 024(a-l  .905)-2540(a-l  .905) 

In  the  above  equation,  it  is  observed  that  the  viscosity  is  varying 
with  the  radius  of  the  indentcr  which  does  not  agree  with  the  theoretical 
results. 
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Moreover  it  is  noticed  from  the  above  equation  that  the  results  will 
improve  with  highly  viscous  fluids  and  this  will  agree  with  the  assumption 
in  deriving  the  theoratical  result,  that  the  viscous  effects  are  dominant. 
Further  in  deriving  the  theory,  the  indentation  of  an  infinite  half-space 
was  considered  while  the  experiments  were  conducted  using  very  finite 
containers  thereby  introducing  end  effects. 

Although  efforts  were  made  to  keep  the  temperature  of  the  room 
constant,  the  experiments  should  have  been  conducted  in  a  precisely 
controlled  environment  to  reduce  the  error  due  to  temperature  variations 
which  changes  the  viscosity  of  the  fluids  to  a  great  extent.  Though  the 
travel  of  indentation  was  fixed  with  the  help  of  supports,  even  then 
variations  were  observed  with  different  weights  of  the  indentor. 

Further,  the  author  performed  the  experiments  to  find  the  viscosity 
of  the  same  test  samples  using  Brookfield  Synchro-Lectric  Viscometer.   It 
was  found  that  the  results  were  quite  satisfactory  within  reasonable  limits 
of  accuracy.  But  the  viscometer  could  not  be  used  at  higher  speeds  with 
viscous  fluids  and  is  limited  for  use  with  fluids  having  viscosity  up  to 
10  poises. 

Recomendations  For  Future  Work: 

It  is  felt  that  if  fluids  of  higher  viscosity  are  used  for  the  purpose 
of  experiment  there  is  every  likelyhood  that  the  results  would  be  within 
very  reasonable  limits  of  accuracy.   It  is  further  observed  that  if  larger 
containers  are  used,  error  will  be  minimized  to  a  great  extent. 
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APPENDIX  A 
FORTRAN  Program  for  Evaluation  of  Co-efficients  by_  THE  METHOD  OF 


LEAST  SQUARES. 
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Notation: 

LI  Fluid  (1) 

12  Fluid  (2) 

13  Fluid  (3) 

R  Radius  Of  Indentor 

W  Weight  Of  Indentor 

V  Velocity  Of  Indentation 

RW  Radius  x  Weight 

RV  Radius  x  Velocity 

VW  Velocity  x  Weight 
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The  mechanics  of  fluids  is  dependent  upon  their  fundamental  property 
named  as  viscosity.  The  exact  information  about  viscosity  is  essential  in 
design  of  pipelines,  in  injection  moulding,  extrusion,  spraying,  bearings 
etc.   Industries  which  deal  with  processing  of  fluids  and  their  by  products 
have  often  faced  difficulties  in  design  of  processing  system  and  handling 
equipment  as  information  available  about  viscosity  of  highly  viscous 
fluids  is  inadequate.  Though  the  commercial  viscometers  are  available  for 
determining  viscosity  but  their  ranges  are  very  limited  and  can  not  be 
used  for  highly  viscous  fluids.   In  light  of  this  a  viscometer  has  been 
designed  end  fabricated  in  Kansas  State  University  which  is  based  on  the 
theoratical  analysis  given  by  Goudy  and  Kirmser.  The  viscometer  is 
designed  for  testing  the  viscosity  of  highly  viscous  fluids.  To  varify 
the  satisfactory  working  of  this  viscometer  the  author  performed 
experiments  for  determining  viscosity  of  fluid  samples  which  were  obtained 
from  Brookfield  Engineering  Laboratories.  The  experimental  data  was 
analysed  statistically  on  computer  and  results  obtained  were  compared  with 
those  obtained  by  Brookfield  Viscometer.  The  results  of  two  viscometers 
show  large  deviation.  The  apparent  discrepancies  are  mainly  attributed  to 
the  small  size  of  the  containers  which  directly  contradict  the  assumption 
made  in  development  of  theory  based  on  indentation  of  semi-infinite  fluid 

surface. 

Results  of  this  study  indicate  that  if  proper  size  of  container  is  used 
the  results  of  two  will  be  quite  close.   The  small  containers  introduce  error 
due  to  change  of  end  effects  which  are  mainly  due  to  propagation  of  the 
disturbance  to  the  container  walls.  This  could  be  avoided  if  highly 
viscous  fluids  are  used  for  determining  the  viscosity.   It  is  hoped  that 


viscometer  designed  vrill  give  better  results  even  for  highly  viscous  fluids 
where  Brookfield  Viscometer  can  not  be  used. 


